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ABSTRACT 


The inherent «elf aligning capability and superior resistance to thermal 
distortion, of the gas-lubricated til ting-pad thrust barring, could often offset load 
capacity limitations and make it an attractive compromise for many applications. A 
numerical treatment for analyzing die til ting-pad thrust bearing is presented along 
with extensive perf or mance information for a particular pad geometry. An optimum 
hydrodynamic crown profile and related dimensions are determined. The range of 
validity of incompressible theory and the rectangular approximation are also 
developed. 
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x x - pad irnkJe radius, in. 
r m * mean radius, in. 
r 0 * pad outside radius, in. 

Tpy * radius to pivot (center of ptwwre), in. 
r ' ■ radius to point where pad is tilted, in. 

T * dimensionless time * — j- 

t * time, seconds 

W * pad load. lbs. 

W t ■ total bearing load. lbs. 

a 3 tilt about radial line, rad 

ag = tilt about axis normal to radius through point of tilt. rad. 
a r 3 tilt about radial line. 3 a? ad 

7 * angle between leading edge and pivot (center of ptr^ni;\^id. 

6 3 crown height, in. 

f 3 radius of curvature of pad. in. 

9 3 angular coordinate, rad 

dp 3 angular extent of pad, rad. 

9* 3 angle to point where pad is tilted, rad 

0 O 3 angle to leading edge of pad. rad 

r 0 2 

A 3 compressibility parameter 3 7 — 

p 4 ° 

H 3 absolute viscosity, Ib«sec/in2 

w 3 shaft speed, rad/sec 

<*> | 3 reference speed 3 ^ , rad/sec 
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1. INTRODUCTION 

Dishing or crowning of bearing surfaces due to viscous heat acneration can 
seriously degrade load capacity of single surface thrust bearings. ^ In addition, 
complex alignment schemes (gimbals, etc.) are often necessary to compensate for 
the distortion of the housing to which the thrust bearing is attached. The 
self-equalized tilting pad thrust bearing with inherent self-alignment and reduced 
tendency to thermally distort, then becomes an attractive compromise. Use of this 
bearing with compressible fluids has been retarded, ostensibly because of poor load 
capacity, but when comparisons are made with imperfect geometries resulting from 
thermal or structural deformations, the tilting pad bearing could be a better choice. 
Also, comparisons are generally made on the basis of equal minimum film thickness, 
an unfair criterion to the tilting p-d bearing, since automatic adjustment of the tilt is 
an added safety factor; i.e., safe operation can be achieved at lower minimum film 
thickness than with the single surface variety. 

This paper presents the pertinent results of investigations of the sectored 
crowned profile pad. Included are optimization criteria, design information and 
procedures and limitations of some of the common approximations. 


Superscripted numerals identify references. 
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2. THEORETICAL FOUNDATIONS 


For isothermal films the compressible Reynolds’ equation in polar coordinates 
is 

a f rph 3 dp | . 13 | ph 3 3p | _ fi a(ph) 3(ph) 

37 37 j + 7 37 3ej 6 se + 12 3t (1) 

By a change of independent variable and non dimensionalization the following 
equation results: 



Q > p2 has been selected as the independent variable because for configurations 
such as a Rayleigh Step, that contains discontinuities in the film thickness 
distribution, it is a more convenient independent variable than the usual P 2 H 2 . 
Although Rayleigh Steps are not considered here the desire for generality of the 
final computer program dictated the use of P 2 as the independent variable. 

Paradoxically, the time dependent terms are maintained in equation (2), 
though the immediate interest is for steady-state information. The reasoning is that 
regardless of whether the time dependent terms are maintained, the resulting 
equation is non-linear and some sort of iterative scheme is required to obtain a 
solution. Steady-state data is derived by diffusing the pressure field as a function of 
time until variations in pressures between successive time intervals do not differ by 
more than a specified truncation error. The advantage of preserving the transient 
terms in the governing equations is that steady-state information is generated at no 
loss in efficiency and the option of conducting future dynamic studies is available. 

Another significant item concerning equation (2) is the use of the variable 
«•>/<•> ]- The compressibility parameter A and time T are non-dimensionalized with 
the reference speed wj, rather than the slider velocity u. If u was used for 
non-dimensionalizing, then ail terms of the right hand side of equation (2) would be 
obliterated at zero slider velocity, and it would not be possible to obtain pure 
squeeze film effects. Thus, with slider velocity, cj/cjj equals unity, otherwise the 
parameter equals zero. 
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Equation (2) is expanded by finite difference and combined into a | 
matrix equation of the form( 2 ) 

tCj] {*,) + IDjI {Q^} + (EjJ <Q J+1 > . { Fj ) 

J - 1, * 

For interior points the [Cj] matrix is tri-diagonal and [Dj] and [Ej] are 
diagonal. The non-linearities occur because the coefficient matrices and the right 
hand side vector are functions of the J~Q. 

The boundary conditions are satisfied by manipulation of the coefficients and 
the right hand side of equation (3). Thus, for ambient pressures on the boundary 

[D] - [E] - 0 
(F> - {1} 


The coefficient matrices of equation (3) are also functions of the clearance 
distributions which depends upon the type of convex profile. Section 5 describes the 
three hydrodynamic crown profiles examined. With regard to film thickness 
distribution, the following expressions apply (See Figure 1 ). 



z 8 + H ‘ + -2 r [ R i coe ( 8 ' - e J - R j - -r* R i * ln ( 9 ’ * «i) 

■HH' 


Y 




Fig. 1 Geometric Variables for Generation of Film Thickness 
Distribution 
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U -».,>>)(?)’ *(«•)’ 

c/c _ ? u 

where 

R m - R m /<CO>9 p/2 ) <5-d) 


H' - c/c 



R' 2 + 2R.R' cos (6 
m l in j 


6 p/2 



For the cylindrical and spherical profiles the radius of curvature enters into the 
Film thickness expression, while for the parabolic it does not appear. This is manifest 
in the final expression for film thickness by the variable Tq/c. Fortunately, as 
demonstrated later, the results are extremely insensitive to the effects of this 
variable so that it does not complicate parametric studies. 


c is defined as the original clearance above the highest point or line of the 
crown, with the pad in the uninclined position. It is used as a 
non-dimensionalization parameter. For any particular set of operating conditions it 
is convenient to maintain c a constant so that non-dimensionalization parameters 
are unaffected by variations in operating film thickness. The parameter Z$ is used to 
vary the uninclined film thickness without changing c. 


The column-matrix method of solution of equation (3) has been treated in 
reference 2. It basically consists of letting 


{ Vl } " lA j ] {Q j } + {B j } 

j - N-2 

Substituting (6) into (3) and solving for iQj <• produces 

( Qj 1 - - [tCj] + t»j1 [Ajlj-Vjl <Q J+1 > 

- - IT | -1 [[Djl {Bj } - tFjlj 

where 

[T] - [Cj] + IOj] [Aj] 

Comparing (6) and (7) 


tA J+l ] 

- - [T]" 1 

‘V 

'V 

- - IT]’ 1 

h ] ,b 3 


( 6 ) 


(7) 


( 8 ) 
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Note from boundary conditions (Q 1 > - {1} which implies that 


[Aj] - 0, (B 2 ) - {1} (9) 

Equation (9) is the starting point for the recursion relations of equation (8), and 
subsequent formation of the matrices [A] and [B]. 

For each time step the non-linearity introduced by the parameter/ q in the 
coefficient matrices and right side is removed by computing its value' from the 
previous iteration, and then solving the resultant linear equation. A simplified flow 
chart for the solution process is shown on Figure 2. Equations (10) through (1 3) are 
for thrust load capacity, righting moments, and viscous friction drag respectively. 



(P-1) Rd6dR 


( 10 ) 


M 

x 



(P-1) R 2 sin 0d0dR 


(11) 


H - f r (P-l) R 2 cos 6d0dR 

y J *V° 


The effect of the viscous drag on the shaft is evaluated from 

* - \ ffft *de<» + ^- f f Y # de<ut 

T i/ T o J ° 1 j V r o J ° 


( 12 ) 


(13) 


The above equations are for one pad. When multiple pads are considered each 
pad is handled separately and the results consecutively added. 


I the franklin institute research laboratories 


SIMPLIFIED FLOWCHART 



Fig. 2 Simplified Flow-Chart Computer Program 
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3. GEOMETRIC VARIABLES 

Figure 3 shows the pertinent geometric variables concerned with the 
subsequent information to be presented. The parameters Tpy and 7 locate the 
position of the pad pivot point. For any particular operating condition this point 
corresponds to the center of pressure. Note that this point is different from the 
origin of tilt (r* .0* ) which is prespecified and required to aenerate the inclined film 
thickness distribution. The center of pressure is necessarily obtained after 
determination of the pressure field and is thus program output. For all future studies 
the pad Was inclined about a radial line through the center of the pad by an amount 
a ■ o r . (tilts in ag direction were not studied. See Figure 1). Also, the reference 
clearance c is redefined as the film thickness along the mid radius at the specified r' 
prior to inclining the pad. (The variable Z, (Figure 1) is now absorbed in c.) 

Some dimensionless parameters are defined in two ways. When investigating 
optimum crown height then it is convenient to maintain c a constant and vary 5 . 
Then c is used as a reference length. When investigating performance of a specific 
pad geometry over a range of film thicknesses., then the crown height 5 is used as a 

ar . 

reference length and replaces c in such parameters as A , — jj— and — etc . 
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4. ANALYTICAL VALIDATION 


The analysis wet checked against published information, and against results of 
an independent analysis using different techniques. Gross^ 3 ) produced results for a 
rectangular slider with a convex profile. To approximate this situation with the 
present analysis, that considers a pie shaped sector, it was necessary to use large 
inside and outside radii and maintain the difference between them at the correct 
slider width. The tabulation shown on Figure 4 shows the comparative results of the 
pressure in the x direction at v * 0. Also shown is the comparative results of total 
load. The comparisons are excellent. As a sector pad approaches a rectangle the Him 
thickness distributions for the parabolic and cylindrical profile should approach each 
other. Indeed, computer runs were made for both these profiles with identical 
results. 

The second check was made against another numerical analysis that computes 
performance of rectangular sliders. Identical grid sizes were used for both cases. The 
numerical approach is different however, in that the rectangular pad method iterates 
on an explicit equation for the pressure distribution. The dashed lines on Figure S 
are pressure variations for the explicit solution, while the solid lines are results from 
the matrix column (implicit) method. Again for the present analysis that treats a 
sector, the rectangular shape was obtained by using radii. Along the row I * 8, 
the pressure for both analy*s folio* identical distributions except at the peak 
pressure points. Along the column J - 8. which is the mid column, the pressures are 
identical. At the column J ■ 14. the explicit analysis gives slightly lower pressures. 
The explicit solution converges asymptotically, while the matrix column approach 
straddles the solution and converges much more rapidly to the coned solution. 
Irrespective, the comparisons are good and the correctness of the program sub- 
stantiated. 
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CHECK CASE (FLAT SQUARE SLIDER) 
COMPARISON OF PRESSURE PROFILES AS A 
FUNCTION OF GRID COORDINATES (I,J) 



Comparative Results, Explicit Solution of a Flat 
Slider Versus Implicit (Matrix Column) Method 
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5. HYDRODYNAMIC PROFILES 

To provide a converging wedge at startup, and to improve load capacity it it 
advantageous to shape the sector with a convex crown. 

Three crown profiles were examined: 

(a) Parabolic - The crown follows a parabolic curve along a circumferential 
line. The maximum crown height occurs along a radial line through the 
center of the pad and is the reference crown height for this configuration. 
When the pad is at zero inclination angle the crown height along any radial 
line does not vary. (Figure 6a) 

(b) Cylindrical - The cylindrical crown is formed by removing a pie shaped 
sector from a cylinder. The crown height increases along s radial line from 
inside to outside radius. The height at the geometric center of the pad is 
taken as the reference crown height (Figure 6b). 

(c) Spherical - The spherical profile is formed such that the four pad comers 
lie on the same latitude circle. The pole position is taken as the location 
for the reference crown height (Figure 6c). 

Film thickness expressions for these crown profiles were described in equation 5. 

Studies of the profiles were made for the two conditions of primary interest 

(1) Performance variations for constant crown height 5 and varying film 
thickness, c. 

(2) Performance variations for constant film thickness ,c and varying crown 
height, 5 

Figure 7 shows load, friction, and minimum film thickness as a function of the 
pivot position, for condition ( 1 ). Note that A is defined with the constant parameter 
6 which replaces c for non-dimensionalization. Three sets of curves are shown tor 
varying values of 5/c. which indicates varying film thickness c. As5/c is reduced or 
as c increases the load capacity decreases, the minimum film increases, and the 
viscous friction decreases. These results are as expected. The cylindrical profile has 
superior load capacity and film thickness. Figure 8 shows Jesuits of condition (2). 
Since c is a constant, A as defined with 6 varies as 6 /c varies. Again, the cylindrical 
profile shows superior load capacity and film thinner. Figure 9 is an expanded 
plot of a set of curves on Figure 8. It is typical for both conditions (Figures 7 A 
8) and it demonstrates conclusively the advantage of the cylindrical profile. For 
equal load coefficients the cylindrical profile provides a thicker film than either the 
parabolic or spherical crowns. The spherical crown exhibits slightly less viscous 
friction hut this is insufficient to challenge the advantages of the cylindrical profile. 
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CROWN PROFILES 



(a ) PARABOLIC CROWN 
8—| 



(b) CYLINDRICAL CROWN 



(C) SPHERICAL CROWN 


Fig. 6 Crown Profiles 
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TILTING PAD THRUST BEARING - PERFORMANCE CURVES FOR DIFFERENT 
CROWN PROFILES, CONSTANT CROWN HEIGHT, VARYING FILM THICKNESS 
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Fig. 7 Crown Profile Performance, Constant Crown Height, 
6 Varying Film Thickness, c 
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Fig. 8 Crown Profile Performance, Constant Reference 
Clearance c. Varying Crown Height, 6 
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U. COMPARISON OF DIFFERENT CROWN PROFILES 

11 PERFORMANCE PARAMETERS vs PIVOT POSITION 




PIVOT POSITION, y/dp 


Fig. 9 Comparison of Different Crown Profiles, Performance 
Parameters vs. Pivot Position 
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Figure 10 are isobar plots of the three crowns at a particular value of A and 
inclination. The parabolic crown has the smallest film thickness, while the 
cylindrical has the largest. The ratio of load to minimum film however, is much 
larger for the parabolic than for the cylindrical profile. Because of the demonstrated 
advantages of the cylindrical profile, further studies were restricted to the cylindrical 
crown profile, except for comparisons between cotpressible and incompressible 
theory, where it was more convenient to use the parabolic profile. 


- 18 - 


Hl^HE FRANKLIN INSTITUTE RESEARCH LABORATORIES 


PRESSURE DISTRIBUTION COMPARISON 
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Fig. 10 Isobar Plots of Crown Profiles at Equal 
Inclinations 
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«. VARIATION OF PERFORMANCE WITH PARAMETER tJ c 

The film thickness expression for the cylindrical and spherical crown profiles 
contain the parameter rjc or rjb if the crown height is selected as the reference 
clearance. The parabolic profile does no: contain this parameter. Examination of the 
comparative results of the three profiles however, do not reveal a marked difference 
among them. The significance of r^c then becomes suspect, and if it can he 
removed as a variable, parametric studies become less complicated 

Studies to determine the effect of Tq/S were conducted for various values of A 
crown height and pad tilt, and Table I is a sampling of results. It ts evident that 
performance is very insensitive to the variable Tq/S. and it need not l*e considered 
significant for parametric studies. 
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7. POSITION OF PIVOT 

A by product of the pie sector pad analysis is determination of the radial 
position of the pivot. The rectangular approximation eliminates this parameter and 
the location is :«suaUy taken at the mean radius. The results of the analysis prove 
that the radial position is very insensitive to film thickness, crown height, inclination 
and compressibility parameter A. Table 2 shows some tabulated results. The largest 
variation is 5%, and the mean radial position (R ro » 0.75) is not very far from the 
correct value. The actual position is slightly outboard of the mean radius. Although 
the hydrodynamic velocities are greater toward the outer radius, the restriction to 
flow is less because of the greater escape area. The counter-balancing effects 
maintains the radial position of the pivot near the mean radius and limits the 
variation to a narrow range. The radial pivot location should be satisfied for 
normal operating conditions so that a pure pitching of the pad is effected. Because 
of the insensitivity of the radial pivot position, rolling of the pad at off design 
conditions should nor be significant. 
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8, OPTIMUM CROWN HEIGHT 

A primary geometric parameter needed to acted a pad configuration is the 
value of the crowd height A senes cf ram were made at different values of the 
co m pr es sibility parameter A maintaining the reference clearance c a constant and 
varying the crown height 5. Since c is the constant parameter it is used in the 
definition of A, rather than 6 . The values of A investigated were 25,100 and 250. 
and 6/c was varied between 0 and 2. This repr e s ents a broad range of parameters and 
sufficient information is obtained to make decisions regarding optimization of the 
crown height Figures 11, 12, and 13 show the variation of Load Coefficient 
Friction Moment Factor and Circumferential Pivot Position as a function of 6/ c. The 
Minimum Film Thickness is used as a family parameter. Since the Minimum Film 
Thickness is p ro gr am output the information in Figures 11.12, and 13 could not be 
obtained directly but required a considerable amount of cross plotting. The 
significant conclusions of these results are as follows: 

(1) Except at very small minimum films, performance is rather insensitive to 
crown height for values of */c between 0.8 and 1.5. 


(2) The carcum: trential pivot position is almost invariant with the value of the 
crown height 

(3) A good selection for the value of tic would be unity, since it provides 
near optimum load capacity over a broad range of film thicknesses and 
values of the compiesnbiLty parameter A. 

(4) As A increases the pivot position has a tendency to move doser to the 
trailing edge and the range of values of the pivot position is reduced. Pivot 
positions between 0.5 and 0.6 give reasonable values of minimum fUr 
thickness. 


(5> A crown considerably improves performance over t flat slider as can be 
seen by noticing (he difference in performance between 5/c ■ 0 and 
4/c ■ 1. For example, for H ^n * .6, the load capacities are at least 
doubled, and a more noticeable improvement is r-o ted at the higher A 
values. Also, except for very low values of film thickness, the viscous 
friction is slightly teas for the crowned profile than for the flat slider. 
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Fig. 11 Variation of Performance with Crown Height, a • 25 
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9. PERFORMANCE FIELD PLOTS 

Field plot* were developed for the purpose of determining performance of a 
sectored pad with a specific crown height and varying film thickness. The pivot 
position can be pre-specified. or alternatively can be optimized for specific operating 
conditions. Since for these plots the crown height is constant it is taken as the 
reference length in the definitions of A and other dimensionless parameters. The 
plots consist of Load Coefficient C^. Minimum Film Thickness Friction 

Moment Factor. F, and Pivot Film Thickness Hp, as a function of pivot position 
y/# p. The crown height to clearance ration, ife and pad tilt. oXq/S, are family 
parameters. Three values of A were considered 25. 250. and 2500. respectively. The 
field plots provide extensive information for the particular pad geometry considered. 

The curves for A * 25 (Figures 14 and 15) are representative of small A or 
nearly incompressible operation. The operating pivot position is spread over a fairly 
wide range. As the film thickness decreases or Sic increases. variations of load 
coefficient are quite sensitive to pivot position. 

In order for film thicknesses to be maintained at reasonable levels, and insure 
adequate load capacity it would appear that the pivot position should be maintained 
between 0.5 and 0.6. Then, for an operating £ /c * I the minimum film would not 
be less than l A of the crown height Also, the value of 5/c * 1 should be selected for 
maximum load conditions i f variable loads will be encountered. Figure 1 5 shows 
Friction Moment Factor aisd the Pivot Rim Thickness as a function of pivot 
position. The viscous friction also shows a sharp increase at higher values of pivot 
position, due to the reduction in film thickness Pivot film thicknesses are not as 
sensitive as minimum films for any constant value of 8/c, and varying pivot position. 
Again, the more acute variations occur at the higher pivot positions. 

Figures 16 and 17 are the field plots at A « 250. For most applications, this 
value of A is well into the compressible region. The overall field plot is much more 
compressed (narrower range of pivot positions) than for the low A case. This was 
also evident in the varying crown height curves. Figures 11.12 and 13. 

The values of load coefficient are almost an order of magnitude higher than for 
the low A situation reflecting the higher hydrodynamic pressures generated by 
increased velocities. Pivot positions in the range betwe e n 0.5 and 0.6S would appear 
to provide reasonable load capacities and minimum films. The Friction Moment 
Factor is much higher for the higher a case differing by an order of magnitude. 

The field plots for A * 2500 are shown on Figures 18 and 19. These curves 
represent operation considerably into the compressible range to the extent that load 
capacity should not be very much improved over operation at much lower A’s. 
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There is slight improvement in load capacity and minimum film as compared to the 
A = 250 case, but the difference is not nearly as much as that between A « 250 
and A = 25. The improvement is obtained by paying an incommensurate price in 
viscous friction losses which are an order of magnitude higher than for A * 250. 
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Fig. 19 Friction Moment Factor and Pivot Film Thickness vs. 
Pivot Position, Constant Crown Height, Varying 
Film Thickness, A ■ 2500 
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10. VARIATION OF PERFORMANCE WITH 
COMPRESSIBILITY PARAMETER, A 

The variation of load coefficient with .% is shown on Figure 20, at various 
values of minimuni film thickness. The studies were made for a specific value of 
t/c a 1.0. An almost Unear relation persists up to a value of A * 20. then the load 
variation levels off, and begins to asymptotically approach a limiting value. The 
Friction Moment Factor (Figure 21) continues to rise hnearty. however. The net 
result is that at the higher A values small gains in load capacity are obtained at the 
expense of a disproportionate increase in the Friction Moment Factor. 
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Fig. 21 Variatton of Friction Moment Factor with A 
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11. DESIGN PROCEDURE 

Although information has been limited to a specific pad configuration certain 
design procedures become evident In most problems the load, geometry, viscosity 
and speed are specified. 

1 . Compute the Load Coefficient for maximum operating load. 

2. Enter a curve such as Figure 20 and select a A between 1 20 and 250. 

3. Establish H m j n from Figure 20. 

4. Compute c from known value of A. 

5. Use a 5/c value equal to I and thus the crown height 6 equals c for 
maximum loading. 

6. From and c, compute minimum film thickness and insure it is 
r easonable for the application. 

7. Compute Friction Moment Factor from a curve such as Figure 21 . 

8. Esuolish pivot position from field plots or other curves such as Figure 13 
or Figure 14. 

9. Using eld plots determine off design performance by maintaining 
constant pivot position. With various values of S/c or ar^S and load 
determine variation in minimum film thickness, friction moment and 
pivot him thickness. 

1 0. Make adjustments as necessary to satisfy limitations. 

The above procedure is intended as a guide in the use of the design information 
presented. There are many ramifications of the given information which would 
require appropriate changes to the procedures. Also, since some data is presented for 
specific values of A,6/c, ar^c etc., cross plotting may be required or else designs 
limited to the parameters computed and plotted. 

As an example of the above design procedure, consider the following sample prob- 
lem. 


Given 

H = 2" 

r o = 4 ” 
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h 


\ -* 


M - 3.19 X \cr° rcyns { lb ~» e .f. ] -airat2!2°F 
' In. 2 / 

N * 20.000 rpm (w * 2094 rad/sec) 

P a “ *5 psia 
W T » 450 lbs (max) 

Use 6 pads. W * 450/6 * 75 lbs/ pad. 


1. Load Coefficieni C. - j 

P r 
r a o 

2. Select A ■ 250 

3. From Figure 20 H m j n = 0.732 

~T 

* 0.414 xior 3 


0.3125 


5. 6/c * I . therefore 6 * 0.414 x I0* 3 

^min “ ^min x c * - 3 ® 3 x * n * 

7. From Figure 21 at A x 25*' H^,, * 0.732 
The Friction Moment F-*.tor, F x 6.75 

Mf » F x p a r 0 2 6 * .6699 in.-lbs/pad 
Mf (6 pads) * 4.0194 in.4bv 

HJI 

HP - Viscous horsepower = 0QQ = 1.276 

8. Pivot Position 

From Figure 16 at Cl * 0.3125, 6/c = 1.0 

f- - 0.585 

p for 0 - 0.7854 

P 


Y -0.4595 rad - 

Alao -f- - 3.625 

a - 3.625 x — 
r 


26.33* 


0.3747 x 10"' 
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From Table 2. r pv /r 0 -a76?. r^ * 3.0o8 in. 

9. Pivot film thickness 

From Figure 17 at */c» 1 and 2L_ = 0.585 
H p = hp/6* 0.857 P 

h p = 0.3544 x 10- 3 in. 


Design performance for other than maximum load conditions can be obtained by 
entering Figures 16 and 17 and following along a line of constant pivot position. 
Table 3 below summarizes some of the dimensionless parameters extracted from the 
curves and ihe computed dimensional performance. 


TABLE 3 


H min C L 

S/c 

OTg/b 

F 

H P 

*T 

(lbs) 

1.0 

0.250 

0.670 

3.85 

5.62 

1.10 

360 

1.5 

0.155 

0.495 

4.90 

4.00 

1.81 

223 

2.0 

0.130 

0.333 

6.00 

3.50 

2.50 

187 

2.5 

0.070 

0.250 

8.00 

2.22 

3.70 

101 


“f 


C 

Nnin 

(total) 

N> 

oxKT 4 


(mils) 

(mils) 

in-tbs 

(mils) 

(rad) 

HP 

0.617 

0.414 

3.346 

0.455 

3.98 

1.062 

1.197 

0.620 

2.382 

0.748 

5.06 

0.756 

1.242 

0.827 

1.787 

1.034 

6.20 

0.661 

1.654 

1.033 

1.191 

1.530 

8.27 

0.420 


Figure 21a shows the variation of load capacity Wy, Friction Horsepower HP and 
pivot film thickness hp as a function of the minimum film thickness h m j n . 
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Fig. 21a Performance Curves, Sample Problem 
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12. LIMITATIONS OF INCOMPRESSIBLE THEORY 

For small values of A. it is often permissable to use compressible ;.:eocy to 
solve incompressible problems or vice versa. It is useful to know the detree of 
inaccuracy in this procedure, and where the bounds of applicability be. When 
considering the incompressible case, the fluid may or may not cavitate in divergent 
film thickness regions, depending upon the ambient pressure. 

Figure 22 shows comparative results between the compressible and 
incompressible cases for a parabolic profile at two different values of minimum film 
thickness. The incompressible theory, because of the linearity of the •overruns 
equations shows a linear variation with A. At low values of A the compressible load 
coefficient is higher, but at the higher values it falls ofT rapidly compared to the 
incompressible case. Further conclusions that can be drawn are as follows: 

1 . Cavitating incompressible theory is closer to compressible theory within 
the range of A’s acceptable to both theories. This conclusion is somewhat 
anamalous since compressible theory might be expected to be more closely 
associated with the aon-cavitating bearing because film rupture cannot oc- 
cur and pressures can go below ambient However, a cavitating pad with 
liquid lubricant has greater load capacity than the non-cavi taring pad because 
negative pressures reduce lo; j capacity. Thus, the cavitating pad more 
closely resembles the higher capacity compressible pad than the non-cavitat- 
ing case. 

2. The range of useful A’s for interchangeability of theory varies as film thick- 
ness. The lower the film thickness, the lower the permissable value of A. 

3. It appears that values of A less than 50 permit use of either theory. 
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13. THE RECTANGULAR APPROXIMATION 

Mo\t of the published performance information on tttting-pad thrust heannes 
treat rectangular sliders. The transformation from a pie-shaped sector is made by 
assuming the rectangular length equal to the mean circumference of the sector, and 
the width equal to the difference between the inner and outer radii. The validity of 
the approximation was tested for two cases, one incompressible and the other 
compressible. For the incompressible case, separate analyses in the form of digital 
computer progiams exist and precise comparisons were made. A flat pad prrrDe was 
considered. Figure 23 show\ the variation in load coefficient with minimum film 
thickness for the rectangular and pie-shaped sectors. Although A is shown as a 
parameter for the incompressible performance, the linear relationship permits 
elimination of this variable. The load coefficient divided by A is analogous to an 
inverse Somme rfeld Number. The A was retained in this instance to maintain the 
same format as for the compressible fluid pads. A fairly large error is introduced by 
the rectangular sector, which increases as the minimum film Sickness is reduced. 
The error is not only due to the differences in hydrodynamic profile, but also due to 
the fact that for the same tilt the minimum film thickness of the rectangular sector 
is greater than for the pie-shaped sector. Alternatively, at <he same minimum film 
the rectangular sector requires a greater tilt to support an equal load. 

With respect to load capacity, the rectangular approximation is over optimistic and 

thus may be misleading. Closer results can be obtained by computing the minimum 
Him thickness of the actual sector, and then calculating the load on the basis of the 
rectangular information. The circles on Figure 23 represent load vs minimum film 
thickness determined in this manner. The results are reasonably close to the actual 

pie-shaped sector performance, and are slightly pessimistic so that a safety factor is 

included. 

Figure 24 are plots of Friction Moment Factor vs. Minimum Film Thickness for 
the incompressible case. Here the rectangular approximation gives conservative 
results at low values of film thickness, and over optimistic results at the higher film 
thickness magnitudes. Using the pie-shaped sector film in combination with ttic 
rectangular values of friction produces an almost parallel curve to the true sector 
except the predicted friction is too low Figures 25 and 2b show variations in load 
and friction as a function of angul^T tilt. This is another common way of plotting 
results. The pie-shaped and rectangular curves for load and friction are reasonably 
close and on this basis the rectangular approximation b good, but misleading 
because minimum film thicknesses can be considerably different. 

Figures 27 through 30 are similar curves for a parabolic crowned pad operating 
with a compressible lubricant. By using large radii and maintaining the difference 
between the radii equal to the pad width a good approximation to the rcctangulur 
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pad was obtained with the program that analyzes the pie-shaped sector. The section 
dealing with verification of the numerical analysis (Section 4) has already 
demonstrated the effectiveness of this approach. 

The results are similar to the incompressible case. Perh a ps a hig he r value of A. 
would have exhibited more dearly defined variations between trends of the 
incompressible and compressible lubricants. 

The following genera] conclusions can be made concerning the rectangular 
approximation. 

1 . On the basis of actual minimum film thicknesses for the pie-shaped and 
rectangular pads, the rectangular approximation. 

(a) will give over-optimistic results for load capacity. 

(b) will give over-optimistic results for viscous friction at high minimum 
films, and conservative results at lower minimum films. 

2. On the basts of pie-shaped minimum film thickness and rectangular 
performance, load capacity predictions will be conservative and viscous 
friction will be optimistic. 

3. Comparisons using equal angular tilts produce reasonably good results, 
being conservative for load capacity and optimistic for viscous friction. 
These results can be misleading however because minimum films can be 
considerably different 

4. As rjr 0 gets largcT. the rectangular approximation should improve. There 
was no attempt to ascertain these effects, or to determine variations with 
A. and geometry. 
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Fig. 23 Comparison Between Pie-Shaped Sector and Rectangular 
Approximation - Load Coefficient vs. Minimum Film 
Thickness - Incompressible Fluid 
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Fig. 24 Comparison Between Pie-Shaped Sector and Rectangular 
Approximation - Friction Moment Factor vs. Minimum 
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Fig. 26 Comparison Between Pie-Shaped Sector and Rectangular 
Approximation - Friction Moment Factor vs. Tilt- 
Incompressible Fluid 
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Fig. 28 Comparison Between Pie-Shaped Sector and Rectangular 
Approximation - Friction Moment Factor vs. Minimum 
Film Thickness - Compressible Fluid 
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Fig. 29 Comparison Between Pie-Shaped Sector and Rectangular 
Approximation - Load Coefficient vs. Tilt- 
Compressible Fluid 
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14. SUMMARY 

A computerized numerical analysis applied to the tilting pad thrust sector 
operating with compressible lubricants has been developed. Extensive performance 
information relating to a specific geometry has been presented along with 
optimization procedures. Limitations of various approximations have been '•frrififd 
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